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ABSTRACT 


An algorithm pertinent to pattern classification and 
cluster analysis of n-dimensional (n>2) data is presented. 
The algorithm simplifies the dimensionality of the sample set 
by a perpendicular projection of each point onto a plane which 
is determined by T orthonormal vectors. The operator views 
the plane on a computer graphics screen, selects a point on 
the screen with a lightpen, and attempts to achieve separation 
of classes by choosing a ‘direction of movement." Each itera- 
tion, a movement of a single point in a certain direction, is 
actually.a reorientation of the plane into which the data 
points have been projected. These iterations continue until 
the operator decides successful separation has been accomplished. 
The method is neither restricted by statistical distribution 
assumptions nor the dimensions of the sample set. From the 
investigation of seven case studies, this algorithm emerges 
as a tremendous assistance in the analysis of multidimensional 


data. 
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І. INTRODUCTION 


The algorithm presented in this paper is an improvement 
on the separation algorithm presented in the thesis by Gilbert 
Lauzon [1]. He developed a separation algorithm to use in 
pattern classification and cluster analysis of n-dimensional 
data. | 

The algorithm was written in Fortran IV for the XDS 9300 
computer installation at the Naval Postgraduate School. The 
Adage Graphics and Text (AGT) system is linked directly to 
the digital computer. The Fortran callable graphics sub- 
routines in the XDS 9300 allow for communications between the 
two systems. Also, by using the library program "MAD" of the 
AGT system, the lightpen and sixteen function switches permit 
an alteration of the display on the graphics screen. 

The computer program projects each point into a plane 
which decreases the dimensionality of the analysis from the 
original n (»2) to 2. The points on the plane are presnted 
to the operator on a cathode ray tube. The user reorients 
the plane on each iteration by selecting a point and choosing 
a "direction of movement." Using the XDS 9300/AGT-10 system, 
the operator makes the decision whether to proceed with a 
different orientation of the plane or terminate the analysis. 
The human decision is based on the pattern of points presented 
on the screen. The interaction of the man-machine system 


is emphasized. 





Many classification problems can be approached in a 
graphical diagram. A two-dimensional plot of the features is 
drawn. From this, the human must suggest a rule that partitions 
the feature space into regions, one for each category. The 
rule can be a series of straight or curved lines depending on 
the number of categories. Since the decision is based on the 
perception of the operator, the rule is quite arbitrary. 
Therefore, the analyst must insure on a proper number of 
samples in each category. In addition, the human tries to 
determine the rule which according to his subjective judgement 
minimizes the probability of error. In essence, this means 
that the probability of a random observation from a specific 
category falling outside its respective region defined by the 
rule should be made as small as possible. The preceding 
analogy is the basic characteristics that prevail in the 
separation algorithm. The graphical diagram is the display 
screen with the projections from the computations of the digital 
computer. The rule is solely determined by the operator from 
the video presentation of the graphical diagram. 

The first modification incorporated into G. Lauzon's 
algorithm was to insure that the A and B vectors were ortho- 
normal on each iteration through the computer program. 
Basically, this means that the two vectors are perpendicular 
and each vector has a length of one. 

Secondly, the existing "X-Y" video display was changed to 


a more discernible pair of symobls consisting of the "Х-0" 





scheme. This improvement greatly enhanced the ability of the 
operator to visually identify the two groups on the graphics 
screen. 

Thirdly, supplementary video information was added to the 
display on the screen. Each iteration of the algorithm 
involves the selection of a point and a "direction of movement." 
These operations are controlled by a sequence of steps 
utilizing the sixteen function switches and the lightpen. The 
function switches have been numbered left to right and top to 
Bottom for easy reference. A record of the use of these 
function switches is most beneficial to the operator. So, 
the numbers of two function switches have been displayed in the 
upper left hand corner of the screen. The one number records 
the function switch corresponding to the type of symbol 
selected and the second keeps a record of the other function 
switches. The modification was incorporated to merely assist 
the operator and did not contribute to improving the separation 


algorithm. 





II. DISCUSSION OF THE PREVIOUS ALGORITHM 


The previous separation algorithm was presented in the 
thesis by Lt. Gilbert Lauzon [1]. The algorithm is applicable 
to use in pattern classification and cluster analysis of 
n-dimensional data. The computer graphics system projects 
each point of a data set into a olane. The projection is 
not necessarily orthonormal. The method reduces the dimen- 
sionality of the problem from n to 2. 

First, a description of the algorithm will be presented 
for a pattern classification problem. Assume there are two 
classes of objects, Р, and Po. These two groups are comprised 
of L-dimensional vectors of numbers. These numbers may be 
measurements or merely some assigned values of nature. The 
goal of the analysis is to determine the existence of a 
partition of the two groups. Once a separation has been 
achieved and a discriminant function has been determined, it 
can be used to classify new objects by the rule: 

Но Ge 2 is a member of the зе Ра: 


f(Z) <0 > Z is a member of the set Ро. 


The two groups are considered separable if an appropriate 
discriminant function can be specified. Otherwise, the classes 


are inseparable and no discriminant function can be assigned. 





The L-dimensional vectors of the measurements on objects 
from the classes Р, апд Р, will be denoted by X and Y, 
respectively. The algorithm commences by calculating the 
average value, DAV, of each of the L-components of the data 


set. 


| Ը1 22 


DAV(K) = (m1 ( 


Here M and N are the number of vectors in P1 and Раз respec- 
tively. Using the L-component vector DAV, the data points are 


centered as follows: 


X'(J,K) 
Π.Κ) 
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Two arbitrary non-zero vectors, A and B, in Euclidean L-space 
are chosen. А linear transformationis computed for each point 
using the two vectors. This amounts to projecting each vector 


to a point in Euclidean 2-space with coordinates 


мћеге А:2 = 7 akzk， 8:24 = 7 0,2%» Z is an arbitrary point 
К-1 К-1 

from the X or Y set. The coordinates are scaled so the points 

can be displayed on the AGT screen. The procedure determines 


the maximum absolute value of both coordinates as shown: 





TMAX = max (max |А"Х| , пах |А-У| ) 
9 Ї 


- 
శి 
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SMAX = max (max |В.Х| ‚ тах |В-У| ) 
J I 


Each coordinate is divided by the respective maximum value. 
Ta ը (Z) = (A-Z/TMAX, B-Z/SMAX) 


At this point, the video display will appear depicting the 
projections of the X and Y points. The operator selects a 
point and chooses a "direction of movement." For example, 
if the tenth X point is picked and a move right is chosen, 


a new vector, A', is computed. 
СООК) = ACK) + X(10,K)-3 K = 1,...,L 


The new values of the A and B vectors are dependent on the 


"direction of movement" and the selected point as follows: 


Direction of Move New A New B 
RIGHT A+ Z B 
LEFT Α - 7 | В 
ЮР | A B + Z 
DOWN A В - Z 


The computer returns to the linear transformation calculations 


and repeats the procedure using the new vectors. 


The goal of the procedure is a separation of X's and Y's 


into disjoint areas of the plane (screen). The operator 


proceeds using the lightpen and sixteen function switches until 


\ 





an appropriate separation is attained or discontinues the 

analysis. If the two classes have been successfully separated, 
the discriminant function f(Z) can be specified. For instance, 
if the X's and Y's on the screen can be separated by a straight 


line, the discriminant function would be of the form: 
#(2) = а (А•2/ТМАХ) + e (B*-Z/SMAX) + f 


where d, e, and f are scalars. The function need not be 
unique. The pattern classification is completed once the 
discriminant function has been determined. 

The separation algorithm for the cluster analysis problem 
is very much like that in the pattern classification. Normally, 
the set of data is not broken down into two classes. The 
entire set is presented as a set of Y's. The concept of this 
method is the determination of groupings from the sample set. 
It is normally used on data where the goal of the analysis is 
to detect clusters. The final groupings of the data points 
into the appropriate number of disjoint subsets is the decision 


of the operator. Each of the subsets represents a cluster. 





III. DISCUSSION OF THE NEW ALGORITHM 


As mentioned in Section I, the separation algorithm was 
modified using orthonormal vectors, A and B. This modifica- 
tion is important for both the pattern classification and the 
cluster analysis problems. The following paragraphs will 
describe in detail the concept of the orthonormal projection 
algorithm. 

First, let us assume the analysis is conducted on a pattern 
classification data set. The data has been divided into N 
vectors in P1 and M vectors in Po. The Р, set is read by the 
computer as the Y points and the Po set as the X points. The 
data points have been centered by calculating the L-dimensioned 
vector, DAV, by averaging each component of the entire data 
set 

Next, the operator must select the A and B vectors to be 
inserted into the computer memory. If the inner product 
between two vectors is zero, then the two are perpendicular. 
The length of any vector is determined by taking the square 
root of the inner product with itself. For example, the length 


of a vector A, denoted by ||A||, is given by 


MALE = (Α.Α) = ( 
K 


մ ԻԶ Է 


լ 


Now the linear transformation is performed on each point 


as follows: 


աա» 





This transformation is an orthonormal projection of the point 
Z (an arbitrary X or Y) into the plane defined by the ortho- 
normal vectors, A and B. The two sets of projections are 
presented to the operator by the video display on the screen. 
The scaling method described in Section II is performed 
on the data points in order for them to be displayed on the 
Screen. By using only one scale factor, the same perspective 
of the two projections will be retained. This scaling term, 


VMAX, is determined as follows: 
VMAX = max ( max ար: Ն. max (|А-У|,|8-У|)) 
9 I | 


where J = 1,...,M and I = 1,...,N. The final values of the 


coordinates from the linear transformation will be 


Та в (2) = (A+Z/VMAX, B+Z/VMAX) 


The first iteration of the algorithm commences when the 
operator picks a point on the screen and chooses a "direction 
of movement." Actually, the computations do not "move" the 
point selected but reorients the plane upon which the entire 
data set is projected. The reorientation of plane is accom- 
plished by computing two new orthonormal vectors. Suppose 
that the user has selected an arbitrary point Z and desires 
to move it right. The first step in the computation is that 


of recomputing the A vector, 


цагг Ն ОЧИК = ту. Е 


12 





For the projection of Z to shift horizontally to the right, 

the new vector, B', should be the same as B. But this is 
Possible if A' and B' are to be perpendicular. Hence, the 
projection of Z shifts also up and down. We want this vertical 
shift to be minimal. It is best to follow the remainder of 

the concept to finding such a B vector with mathematical 


expressions. 


The appropriate symbols will be defined as: 


A and B = old orthonormal vectors 
A' and B' = new orthonormal vectors 
|А-В| = absolute value of the inner product 


of the vectors A and B 


The new vector A' must be changed so it remains constant with 
a unit length. This is done by dividing the vector (A + Z) 


by its length as such: 
А' = (А + 2)/||(А + 2)|| 


This is the new vector which will be used in recomputation of 
the new projections. 

The new B vector will be determined using Lagrangian 
minimization. This technique is applicable because the require- 
ment to determine the minimum value of the objective function, 
[2:(B' - B)|, has the two restrictions imposed on it. Now, 
to find the B' which minimizes the difference between Z-B' and 


Z*B, the following problem will be solved: 


l= 





Minimize |Iz-(B' - B)| (1) 

subject to В'•А' = 0 (2) 
НЕ 

(| (3) 


Using the Lagrangian method to find the solution, let 
F- € Z-(B'-B) - A4 (B'-A*) - 22 (118115 - 1) (4) 


where Хү апа Ха are the Lagrangian multipliers. Taking the 
first derivative with respect to B' and setting the result 
equal to zero: 


Q E = І і = 
Sei = £ Z - XIA' - 2A?B' = 0 (5) 


Using equations (2) and (5), the value for dy is 
Ay t ZA! (6) 


This result is true if and only if ho has a non-zero value. 
Using equations (3) and (5), the value for А2 сап be determined 


as: 


7 - (7,А') А' 
LAT] (7) 


ДЕ 


2 


Substituting the values of A, and ho into equation (5), the 


final solution to the minimization problem is: 


(eo: «ΑἹ 
С ГОЗ Е ОД P 
The final step is the determination of the sign to use for the 


B' vector. By using the inner product of Z with both sides of 


14 





equation (8) апа the second derivative of equation (4), the 


following statements can be shown to be true: 


D Use + Sign if sign (Z-B") Z sign (7 - (Z-A')) 
(2) Use - sign if sign (Z-B') = sign (Z - (Z-A')) 


The same reasoning can be extended to "moving" an arbitrary 
point Z to the left. However, in selecting a point Z and 
"moving" it up, the new vector B' will be determined first. 


This time the minimization problem will be: 


Minimize zen л 
Subject to В'•А' = 0 
2522 
| |ո'||՛ = 1 


The following table gives a guide to the computations 


of the new vectors using an arbitrary point Z. 


Direction of Move New A New B 
RIGHT (A + Z)/||A + Z]| B* 
LEFT Ա Մայա B* 
UP A* ΟΠΕ ΟΙ 
DOWN A* БЕТІ ||| 


*Denotes the use of the minimization technique. 
These new A and B vectors are used in the recomputation of the 
linear transformation T on the data set. After the projections 
have been calculated, a new display of the plane depicting the 


two subsets will appear on the screen. The operator can make 


15 








another selection or terminate the analysis depending on the 
presentation of the points. 

E "X-Y" letters in the old routine were changed to the 
"X-0" symbols. This was done to make the presentation easier 
for the user to discriminate the two data sets. The similarity 
between an X and Y made it perplexing when the two different 
symbols were in close proximity. The two letters are sketched 
on the screen through a set of draw and move instructions. 
These computer instructions are stored in a vector called IPACK. 

Video information of the function switches was also 
incorporated into the computer program. Now, when the operator 
depresses one of the sixteen function switches, the number of 
the switch is displayed in the upper left corner of the screen. 
The number in the first row indicates the function switch for 
selecting either an X or a Y point. The value of these 
commands will be an 11 or 3, respectively. The remaining 
function switches will be indicated on the second row. This 
additional video display is controlled by a series of ENCODE 
and CALL TEXTO instructions in the computer program. The role 
of the function switches is to communicate between the Fortran 
graphics package "MAD" and the 9300 computer. The display of 
the function switch numbers assist the operator in efficient 


execution of the program. 
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IV. ANALYSIS OF THE SEVEN CASES 


The operation of the separation algorithm using the 
orthonormal vectors was first examined using four sets of 
artificially generated coordinates. The first set was just 
to test the modified program. The remaining three were 
generated for the sake of testing the applicability of the 
method to pattern classification. Each set of points filled 
various geometrical shapes. This quality allowed for ensuring 
proper appraisal of the procedure using multi-dimensional data. 
Finally, three cases of measurements from "real world" applica- 
tions were analyzed to evaluate the performance of the new 


algorithm. The seven cases that will be discussed are: 


TESTS: 





(1) One three-dimensional pyramid 

(2) Two three-dimensional cubes 

(3) Two six-dimensional cubes 

(4) One four-dimensional shell enclosing a four-space 


sphere 


APPLICATIONS: 
(5) Iris flowers 
(6) Fossils 


(7) Coronary heart disease 


17 





Prior to exemplifying on each case, a few comments are 
appropriate. The output from the computer printer comprises 
a sixty-one by thirty-seven dimensional matrix. The printout 
is a six by six inch square. The outline of the square is 
shown by asterisks and the coordinate axes are indicated by 
four zeroes. The symbols denoting the data points on the 
printout correspond to the plot on the AGT graphics screen. 
The points will appear as an X, 0, 8, E, or J. The "X" and 
"0" symbols match the positions of the X and Y projections, 
respectively. ` Ап "8" symbol represents two or more Y's in a 
cell; and "E" refers to two or more X's. A cell denotes a 
matrix element of the printout array. The matrix is stored 
in the computer memory. For example, PIC (31,18) in the 
computer program would refer to the symbol in the cell occupying 
гом 31 and column 18 of this matrix. Lastly, the "J" letter 
indicates a mixture of X's and Y's occupying the same cell. 
The coordinate axes have been drawn in by hand. The encircled 
points signify the X or Y selected in the next iteration and 
the arrow represents the direction of movement. 

For each case, the initial values of the A and B were 
extracted from the table in Appendix IV. This table represents 
only one group of an infinite set of possibilities. This 
particular group was chosen because the absolute values in any 
dimension are nearly equal. The author felt the use of nearly 
equal components would initially bring all measurements into 


the projection by the same relative proportion. It was thought 


18 





that this should give a better intitial presentation of the 
projections than any other choice of values. However, it 
should be pointed out that the decision on the particular set 
of values to insert is left up to the user. 

Polaroid prints of the AGT plot of the projections have 
been included with the computer printout. The close resemblance 
between the printout and the Polaroid pictures is quite 
apparent. Of course, the Polaroid picture is a true repre- 
sentation of the projections. Any discrepancies in the plot 
from the computer printout would be due to the restrictions in 
setting up the matrix array. Adjustments of the vertical and 
horizontal controls were necessary to center the plot on the 
screen so a Polaroid print may be distorted but the relative 
distance between points has been maintained. The elaboration 


and the conclusion of each case will now be presented. 


CASE 1 

The first case involved the coordinates of a three- 
dimensional rectangular pyramid. The points represent the 
outline of this figure. The data set consisted of twenty-one 
vectors. Each consecutive iteration will show the rotation 
of the geometrical shape. The change in the A and B vectors 
merely reoriented the perspective of the figure. Also the 
Shape of the figure remains true with each iteration. This 
example was used in evaluating the successful operation of the 


orthonormal projections. 





Figures IV-1-1 through IV-1-5 in Appendix V show a five- 
step sequence in the application of the orthonormal projection 
concept. Since this case was done only as a test of the 
algorithm, no printout of the appropriate terms is given. The 
outline of this geometrical shape was drawn in by hand in the 
figures. The Polaroid picture of the AGT display, Figure 
IV-1-5P in Appendix V, is included to compare with the print- 


out given in Figure IV-1-5. 


CASE 2 

Two three-dimensional cubes were generated for analysis 
using the pattern classification technique. The X and Y 
subsets comprised twenty-tuples depicting the outline of the 
two geometrical figures. The Y subset has its closest vertex 
to the origin at the point (3,2,5). Each of its sides has a 
length of two units. The other cube (X subset) had its closest 
vertex at (7,6,9) and also two units on a side. The two cubes 
were made disjoint. This made it feasible to specify a dis- 
criminant function depicting the partitioned subsets. Again, 
this was another test in evaluating the procedure but involved 
two disjoint figures. 

Figures IV-2-1 through IV-2-5 in Appendix V show a six step 
sequence in the application of the procedure on the two cubes. 
Separation was maintained throughout the sequence. The values 
of the A and B vectors have been indicated below the figures. 
The initial figure used the values from the table in Appendix 


IV. 
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A possible discriminant function in Figure IV-2-5 can 


be specified as 


f(Z) = A-(Z - DAV)/VMAX 


where A has the values printed from Figure IV-2-5. If 
ШИТИ 0, Z is a member of the subset of Y's; if f(Z) < 0, 
Z is an entity in the X's. These results can be checked by 
the following example. The printout indicated the values 


for VMAX and DAV; 


DAV(1) = 6.000 DAV(2) = 5.000 


DAV(3) 8.000 VMAX 4.333 


Assume that Z has the coordinates of the point Y(1) in 
Appendix III. The calculations indicate the number for f(Z) 
is + 0.662. The computation has shown that Z is, in fact, a 
member of the subset of Y's. Also any other point from the 
Y's will satisfy the condition that f(Z) will be positive. 
One of the two cubes in Figures IV-2-4 and IV-2-5 have 
been drawn in by hand. Also a figure of the Polaroid print 


is included to compare with Figure IV-2-5. 


CASE 3 

For the third application using the new separation algo- 
rithm, the dimensional space of the cube was increased from 
three to six. Seventy-tuples were generated for each subset, 
X and Y. The unit cubes were generated on the IBM 360 computer 


as uniform distributed random numbers. A point from the Y 
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subset would be denoted by Y(I) 1-2) Y AS), 
Y(1,4), ИН. (I 6)) i wheremi =o. N and 1 < Y(I,K) 
ա - 1,...,6; and for the X subset X( J) = (X(J;,1]), 20 
MS, XxX(J,4), X(J,5), X(J,6)) where J з 1,...,М and O < 
աի 68, Է-|....,6. Again, the two cubes were generated 
so the two samples were disjoint. First, if one of the itera- 
tions showed a combination of "X's" and "O's" in an area of 
the plane, it would be easier to separate the two groups on 
subsequent steps. Secondly, once separation between the two 
groups was attained, a discriminant function could be defined. 
The printout of the graphical display is shown in Figures 
IV-3-1 through IV-3-6 in Appendix V. The degree of complexity 
for the two samples was increased from that of Case 2, but 
a partition of the two groups was still achieved. Since this 
case is so similar to Case 2, the discriminant function will 
not be determined. Consequently, the A and B vectors will 
also be omitted. In Figure IV-3-2, the two cubes shown by 
their projections can not be separated because of the mixture 
of "X's" and "O's" around the center of the graph. It was 
easy, however, to separate the two as shown in the next step, 
Figure IV-3-3. The similarity between the printout and the 


Polaroid snapshot is shown in Figures IV-3-5 and IV-3-5P. 


CASE 4 
The previous two cases were applications of the parti- 
tioning of geometrical shapes which were separable by a hyperplane. 


A sphere and a shell in four-space were created for the purpose 
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of making the division of two groups more complex. Twenty- 
tuples were generated for the X subset and forty-tuples for 
the Y subset. A typical position of a vector in the Y set is 
ПОНЕ (УСТ). 111.2), Y(1,3), Y(1,4)), where I = 1,...,М. 


The components must satisfy the requirement 2 га + Y(1,2)* 


2 + Y(1,4)? « 4.. Similarly, an X point must fulfill 


3) 
the condition that X(J,1)% + x(9,2) + x(9,3)% + x(2,4)* « 1. 
The two subsets are again disjoint and they form a sphere 
surrounded by a shell separated by one unit in radial distance. 
Two different trials on the sample set were conducted. The 
first run was finally terminated after twenty iterations without 
achieving successful separation. It should be noted that the 
arrangement of the two groups made it impossible for an appro- 
priate partition. Of course, this is assuming that the shell 
has the proper density of. points and no voids exist in it. 
However, if the procedure is modified, separation can be 
feasible. The second run was carried out by including the 
Square of each component in each entity. The dimension of the 
Sample set becomes twice the original. To clarify this last 
remark, if there was an arbitrary point, Z = [ప іп 


the original set, it will now become Z = DELI ԱՆՆ. 


B 2,5. The second trial was conducted using this eight 


component vector on both sets. 


Z 


The first trial on case four appears in Figures IV-4-1 
through IV-4-5 in Appendix V. These figures confirm that the 
separation was never attained. Only the first six steps have 


been included in the figures. The remaining steps were omitted 
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since they did not improve the pattern of the projections. 
Thus, a discriminant function can not be specified.  Conse- 
quently, the printout of the A and B vectors will be eliminated 
from the discussion. 

The situation improves in the second run using the eight 
component tuples. This is depicted in Figures IV-4-6 through 
IV-4-10 in Appendix V. Separation was attained in a mere three 
Steps of the algorithm. This modification of the procedure 
can be attempted on samples that can not be successfully 
partitioned using the basic method. The one subset may be 
enclosed inside a large sphere-like shape including the other 
Subset. This is one approach to obtaining a division of the 
two groups. The components of the A and B vectors in the 


second run are printed on the bottom of each of the figures. 


CASE 5 

This sample set represents the first of three "real world" 
cases. The data set for this case portrays multiple measure- 
ments on three species of Iris flowers. The data was originally 
used by Ronald A. Fisher in his classic paper on discriminant 
analysis [3]. It has become widely used in illustrating cluster 
analysis procedures. The four measurements collected are 
(1) sepal length, (2) sepal width, (3) petal length, and (4) 
petal width. The sample set has been reduced from the original 


150 to 75 samples for this particular analysis. 
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This particular sample set has the characteristics for 
either a cluster analysis or pattern classification problem. 
The author will approach the investigation using the pattern 
classification technique. It is quite possible from looking 
at the values of the three species that the one, Iris setosa, 
will have projections in an area of the plane by itself. The 
projections of the other two species appear to be nearly the 
Same. Hence, it may be somewhat сше to partition these 
latter two species. Therefore, the observations from Iris 
versicolor and Iris virginica will bedivided into two different 
subsets. Тһе Y subset will be composed of Iris setosa and 
Iris versicolor; while the X subset of Iris virginica. Each 
species of flowers was composed of twenty-five vectors in 
four dimensional space. 

The complexities in measurements is apparent in the parti- 
tioning of these three groups. Outliers may exist which 
indicate an observation is far removed from the main body of 
the data. They result from the fact that the sample under 
investigation may contain a certain proportion of extreme 
cases or from a population other than the sample from which 
the rest of the data came. An observation from the Iris 
versicolor species presents itself as a possible outlier in 
@aemanalysis of the graphic plot. This point is Y(28) in 
Appendix III which coincides with the projection of X(24). 

The two appear as one point in the last two figures as a "J" 
symobl. The point Y(28) could not be separated from the group 
of X's as it is projected to the left of the line denoting 


the discriminant function. 
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Figures IV-5-1 through IV-5-6 in Appendix V show the 
twelve step sequence in the analysis of the three Iris 
flowers. Even though a discriminant function is generally 
determined in partitioning two groups, two linear functions 
will be specified which divide the three types of flowers. 
This will be done using Figure IV-5-6. The first function 
will separate the. X subset from the Y subset. Then the 
second will divide the two Iris species comprising the Y 
subset. First, the values of the A and B vectors from the 
bottom of the figure must be recorded. 

Now, a linear discriminant function has the form given in 


Section II as; 


e(z) = SA(Z - DAV)) , 


VMAX VMAX 


where d, e, and f are scalars. The f scalar is the negative 
sign of the ordinate intercept while the term d has the 
negative sign of the slope of the line. Finally, the value 
of the coefficient e will be one. Using the line in Figure 
IV-5-6 separating the X and Y groups, an appropriate dis- 


criminant function has the form: 


_ 6(A(Z - DAV)) , B(Z - DAV) 
(2) VMAX t mag AL ls 


The magnitude for DAV and VMAX from this figure are: 


DAV(1) = 5.872 DAV(2) = 3.061 
DAV(3) = 3.804 DAV(4) = 1.212 
VMAX SUIS 
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Using these figures, the value for f(Z) is - 0.054 when the 
components of X(24) from Appendix III are substituted into Z. 
Also carrying out the same computations, the result for f(Z) 
is * 0.252 for the point Y(26). These calculations have 
merely emphasized that a positive value for f(Z) implies Z 
is a member of the Y subset. If f(Z) is negative, then Z is 
an entity of the X subset. 

Now, a description of the line splitting the Y subset 
must be resolved. The form of the discriminant function has 


been determined as: 


f(U) = A UDA 2 0.370 


VMAX 


where the number 0.370 is the value of the intercept on the 
abscissa axis. If the coordinates of Y(38) are substituted 
for the arbitrary point U, the result of f(U) is - 0.030. 
This calculation indicates that any point from the set of 
Y's, where f(U) < 0, represents an entity from the Iris 
versicolor species. Any U where f(U) > 0 corresponds to an 


observation from the Iris setosa. 


CASE 6 

The sample set for this analysis problem is quantitative 
measurements of morphologic parameters of specific fossils. 
The data set was extracted from Wright and Switzer's paper 
titled, "Numerical Classification Applied to Certain Jamaican 


Eocene Nummulitids" [4]. Eighty-two eight component vectors 
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were used in the investigation. Six of the original observa- 
tions had missing parameters and were therefore omitted from 
ене sample set. 

Again like Case 5, this sample will be divided into two 
groups for a pattern classification problem. The author feels 
the separation will be easier to attain using two subsets. 

It should be noted that Lauzon investigated these specimens 
using cluster analysis procedures [1]. The Y subset was 
initially composed of the first thirty-eight observations. 
It is comprised of specimens from the Operculinoides wilcoxi 
species. This species will be denoted by Group A. The X 
subset was composed of the other three species: Operculinoides 
jennyi, Operculinoides cushmani, and Nummulites striatoreti- 
culatus. These three types of Nummulitids will be called 
Groups B, C, and D, respectively, Using these subsets, it 
was extremely difficult to pick out the projections from 
Group B and those from Group C. In order to eliminate this 
confusion, the two subsets were rearranged. Observations 
from Groups A and B were combined to form the Y subset; 
Groups C and D joined to make up the X subset. The second 
run was conducted using the above groupings. Also to make 
the investigation similar to the discussion in Wright and 
Switzer's paper, only the first six parameters were used. 
The results from the analysis were quite remarkable. Even 
on the first step, three different groups were apparent; 


Group A, Groups B and C, and Group D. In Wright and Switzer's 
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investigation, ninety-six iterations were conducted to deter- 
mine the coefficients of a linear discriminant function. They 
were able to find two different groupings, Groups A and D and 
Groups B and C. The same procedure using 96 iterations was 
repeated to find another linear function dividing Group A and 
D. At gast, a linear function splitting Group B and С was 
determined. Then. the two researchers utilized these two final 
divisions as input to a standard discriminant analysis program. 
The output resulted in two canonical variables. The plot of 
these two variables revealed three groupings: Group A, Groups 

B and C, and Group D. The two concluded that the specimens 
from Groups B and C could not be partitioned. The results 

in this thesis are in agreement with these findings. The 
procedure used by this author appears to be superior in two 
ways. First, there are substantially fewer computations, 

hence the computer time to accomplish the analysis is consider- 
ably less. Secondly, the plot of the results is readily 
available from the printout of the computer. It should also 

be noted that G. Lauzon's investigation of this sample set 
Showed only two groupings. 

Figures IV-6-1 through IV-6-3 in Appendix V show the 
graphical display obtained on the first run. These figures 
depict the eight-dimensional analysis of the fossil measure- 
ments. In Figure IV-6-2, the three groups have been labeled. 
There is no visual method to distinguish between the B and C 


groups. The displays of the second run are given in Figures 
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IV-6-4 through IV-6-7. The values of the A and B components 
of these figures are given below each. 

The two discriminant functions separating the three 
groupings will be specified using Figure IV-6-7. First, the 


tabulation of the numbers for VMAX and DAV are as follows: 


DAV(1) = 163.220 DAV(2) = 41.43 
DAV(3) = 3 ՉԱ DAV(4) = 25.012 
DAV(5) = 88.866 DAV(6) = 677.817 
VMAX Հ 587.954 


A description of the line splitting Groups A and D must be 
determined. The form of the discriminant function has been 


resolved as 


B-(Z - DAV) 


VMAX 202639 


f(z) = 


If f(Z) > 0, then the arbitrary point Z is a member of Group 
А, В, ог С. On the other hand, if f(Z) < 0, then Z is an 
entity of Group D. The equation of the discriminant function 


describing the upper line is given as: 


_  0.6(A+(U - DAV)) , B-(U - DAV) 
f(u) VMAX $ VMAX 


U is an observation from Group A or D if f(U) > 0. Likewise, 


U is a vector from Group B or C if the result of f(U) < 0. 
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CASE 7 

The final analysis was conducted on a thirteen-dimensional 
pattern classification problem. Coronary patients from 
Letterman Hospital comprised the sixty-tuples making up the 
Y subset [1]. The other subset was composed of fifty Army 
personnel at Silas B. Hays Hospital who had undergone annual 
physical examinations. Those patients from Letterman had 
been diagnosed as having some type of coronary heart disease. 
The X subset consisted of thirteen features from a random 
sample of the Army personnel. There were no indications from 
the latter group to either classify or not as having heart 
disease. If successful separation is obtained and a discrimi- 
nant function is computed, then it could be used to classify 
new patients. 

A few of the components in each thirteen-tuple represented 
measurable entities such as weight and age. Other character- 
istics such as race and smoking behavior were assigned 
arbitrary values rather than numbers to signify a measurement. 
The analysis will show that these classifiers are not really 
appropriate to use. The actual measurements exhibit at least 
the characteristics of an interval scale. This is a minimum 
requirement in computing a meaningful average value of a 
parameter. Those components which had an arbitrary value 
assigned do not meet this criteria for an interval scale. In 
fact, they have the characteristics of an ordinal scale. 


Therefore, when the average values of the components, DAV, 
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are computed in the program, it is inappropriate to use these 
classifiers in the analysis. 

The first run was conducted on the thirteen-dimensional 
sample set. The displays of the projections are given in 
Figures IV-7-1 through IV-7-6. The last figure had the best 
presentation with respect to partitioning the two groups. A 
line dividing the-"X's" and "O's" was drawn in the figure. 

If a discriminant function had been specified for this line, 
thirty-six of the fifty X points and forty-nine of the sixty 
Y points would have been properly classified. This result 
means that approximately 25% of the total sample have been 
projected out of their respective regions. This last state- 
ment does not seem satisfactory for any prediction model or 
separation. 

The values of the A and B vectors used by G&G. Lauzon of 
this data indicate a definite pattern. For those assigned 
values, the number in both vectors for each component are less 
than or equal to the absolute value of two. Оп the other 
hand, one of these two numbers is very much greater than 2 
for every measured component. For this reason, the components 
of assigned values will be deleted in the second run. In 
addition, this algorithm is more appropriate for separation 
of data where the sample set is on an interval scale. 

In the second run, an arbitrary observation Z will be 


denoted by 


τι, 12722113) 120154), 200,12), Z(H,13)) 
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where H = 1,...,M for an X and H = 1,...,N for a Y. Figures 
IV-7-7 through IV-7-12 show the printout of the five-dimensional 
plot. The figure with the best partitioning is IV-7-12. This 
time the discriminant function would result in a smaller 
probability of error. There are fifty-four Y projections in 
the left region and thirty-six X's to the right. Now, the 
result is less than 80% of the projections fell into their 
respective areas. The author feels that this second run has 
also led to an unsatisfactory separation of the two groups. 
The reason for this conclusion is the high density of points 
near the dividing line. In addition, the distance is rather 
insignificant for the different points of the two subsets 
across this dividing line. Therefore, the author feels that 
an attempt to specify a discriminant function for this data 


would be futile. 
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у. CONCLUSIONS AND RECOMMENDATIONS 


It can be concluded that the improvements on the algorithm 
have resulted in achieving greater success than the previous 
application. First of all, the revised algorithm provided a 
partitioning of the generated sample sets in the first four 
cases. The new separation algorithm was also applied to three 
"real world" applications. The first investigation obtained 
EDU dactory partitioning of the three Iris species. The 
three groups were divided by specifying two discriminant 
functions. The second sample set involved a more complex 
analysis problem. The observations consisted of an eight- 
dimensional array. In addition, it had been hypothesized 
that there existed four species in the data. The display for 
this analysis indicated the presence of only three species. 
The same conclusion was also reached by Wright and Switzer [4] 
in their approach to this discriminant analysis problem. 
However, their computations in reaching this result were much 
more involved than the algorithm in this paper. The investi- 
gation by G. Lauzon showed the presence of only two such 
separate groupings. From this, the improved algorithm using 
the orthonormal vectors appears to achieve a more definite 
partitioning. Finally, the employment of this method on the 
thirteen-dimensional medical data resulted in unsatisfactory 
separation. The rejection of a partition may be due to the 


use of inappropriate parameters in analyzing coronary disease. 
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All separations are really approximations. The final 
classifications must be subjected to a variety of sample sets 
for a confirmation of the discriminant function. Even after 
this has been accomplished, the presence of outliers will 
result in some finite value for the probability of error. The 
description by the discriminant function is only the optimal 
division of the groupings. 

The use of the values from the table in Appendix IV should 
be mentioned. All investigations in this paper were conducted 
choosing the initial orthonormal vectors as listed in the table 
(Appendix IV). However, if different values had been sub- 
stituted in the components, the final separation may be better. 
In the future, the operator should utilize a variety of unit 
orthonormal vectors in the attempt of achieving a successful 
partition. 

This algorithm has provided a relatively easy and simple 
approach to the classification problem. It is possible that 
the merit of this method could be incorporated into a portable 
unit. A computer module designed for the appropriate computa- 
tions would provide a means for obtaining the required 
projections. A television cathode ray tube could act as the 
display system. An interface system could be developed to 
handle the input/output operations. The outstanding feature 
using a color t.v. set is the ability to display a variety 
of colors. The various subsets of the sample set could be 
linked to a respective hue. The combination of these two 
systems would make up an ideal graphics/computer system in 
which the algorithm would be a major feature. 
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The development of this algorithm could be applied in an 
unlimited variety of fields. The analysis of medical and 
taxonomy samples has already been mentioned. The method is 
pertinent to any field in which there is mutivariate data to 
be separated. For instance, measurements of meteorological 
parameters may be subjected to this method of ANE SSS The 
motivation of the investigation could be to determine the 
existence of conditions pertaining to the development of a 
tropical storm and those present during the infancy of a 
hurricane. This algorithm may be able to distinguish the two 
different groupings. If this can be accomplished, then the 
model can be used as a basis in predicting future storm 
activity. 

The final comments are directed to future changes in the 
computer routine. First, the A and B vectors should be 
punched on computer cards and read in with the data deck. 
This would eliminate the possibility of typing in an incorrect 
component. Secondly, a video presentation of the number of 
the selected point would be most auto the operator. 

In addition, the printout of the plot should indicate the 
selected point by a different set of letters than the present 
"K-0" scheme. These two modifications would make it possible 
to repeat the same iterations on subsequent runs. Lastly, 

it would be beneficial to the operator if a series of sample 
sets could be read. This could be arranged by adding the use 
of another function switch. The operator would now be able to 
conduct a series of different investigations without leaving 


the graphics console. 
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APPENDIX II. PROGRAM OPERATING INSTRUCTIONS 


There exists two methods by which this program can be 
executed on the XDS 9300. The first and fastest method is 
utilizing a magnetic tape. The memory dump of the compila- 
tion from the digital computer has been recorded onto a 
magnetic tape. The elapsed time until execution commences 
would be approximately five minutes using this tape. The 
other method is simply reading the Fortran program from 
the computer cards. It would take roughly thirty minutes 
until the program is ready to execute using the cards. This 
. time factor is very crucial if more than one set of data 
cards must be read. 

In executing the program from the magnetic tape, the 
operator may consult the publication titled "Electrical 
Engineering Computer Laboratory" by R. D. Delaura [2]. A 
brief guide in using the magnetic tape mechanism applicable 
to the object language filed on this, tape follows: 

(1) Load the magnetic tape 
(2) Unit select to position "1" 
(3) Automatic/manual switch to "MANUAL" 


(4) Depress forward drive (tape will stop at 
the LOAD POINT) 


(5) Automatic/manual switch to "AUTOMATIC" 
Ihe magnetic tape is now ready for instructions from the 
Х05 9300. The operator must incorporate five control cards 


into the data deck. These control cards are instructions 
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to the various systems used in the execution of this program. 


The deck of cards will be as follows: 


(1) 8991 
(2) 408 
(3) AGT 


(4) ASSIGN X1 з МТ1А 

(5) RERUN 

(6) (DATA DECK) 
The deck 15 пом ready to load into the card reader and 
execution can commence. 

The following instructions will be applicable to the 

use of computer cards. The list of control cards, that 
must precede the source program printed in Appendix I, 


are given below: 


(1) 8991 
(2) $ DATA 
(3) $PATCH 


(4) 007067 001 06711 


(5) $ЕМО 
(6) 408 
(7) AGT 


(8) FQRTRAN LS, GP 
The data deck is placed at the end of the Fortran program. 
The deck of cards is now ready to load and execution can 
commence. 
In preparing the computer program using either method, 
four cards must be punched and inserted into the data cards. 
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The first three cards can be used for the purpose of verbally 
describing the data to be investigated. This information 

can be punched using any part or all of the eighty column 
card. The operator must insure in insert these three cards even 
though one or more cards may be blank. The fourth card is the 
instruction by which the computer reads the data from the 
punched cards. For illustrative purposes, The statement might 
БЕЛО, " (4F7,2.3X.3F9.2)." This indicates that a total of 
seven parameters will be read from each card. The first 
parenthesis must be punched in column seven of the card. 

The actual data cards are inserted behind these four cards. 
The four cards plus the cards of data will be considered 

the data deck. 

A guide for operating the XDS 9300 computer and the AGT 
console using either a magnetic tape or cards will be given. 
First though, assume that there are M number of X points 
and N number of Y points. The value of M and N must be 
less than or equal to eighty. In addition, all observations 
in both sets constitute an i cimensional wee tor: The number 
of L must be less than or equal to thirteen. Now, the 
instructions will be presented. 

(1) Use the "POWER ON" procedure for the XDS 9300 

and the AGT graphic systems. 
(2) Executing the program 
(a) At the XDS 9300 control panel, depress "SENSE 
2" if using cards. Do not depress if using the 


tape. Then depress "IDLE", "CLEAR" and "CLEAR 
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Е 


(b) 


(c) 


FLAGS" simultaneously, "RESET", "RUN", and 
"CARDS." The card reader will commence reading 
the punched cards. The line printer and 
teletypewriter (TTY) will simultaneously print 
"JOB." The computer will now start compiling 
the program. 

At the time execution of the program starts, the 
TTY will print the message: "TYPE IDEV = 2* AND 
A CARRIAGE RETURN IF USING AGT2." The input 
light will illuminate to alert the operator. 

The operator must type the number of X's and Y's 
plus the dimension of the two subsets of data. 
In addition, the operator has the option of | 
printing all WRITE statements prior to statement 
#57 in the computer program. This is done by 
typing "NPFLG = 1" if a printout is desired; 
NPFLG = 0" for no printout. As an example, 
Suppose there are 15 X's and 20 Y's in 10 
dimensions. If the ne using AGT 2 and 
requests no printout, type "M = 15, N = 20, L = 10, 
NPFLG = 0, IDEV = 2, *" and a carriage return at 
the computer TTY. If using AGT 1, simply delete 
"IDEV = 2." The card reader will begin the 
process of reading the card. 

The word "NAMELIST" and a blinking cursor will 
be displayed on the bottom left corner of the 


AGT screen. This will occur when computer reaches 
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Statement #57 in the program. The operator 
must type in the components of the A and B 
vectors. The values from Appendix IV may be 
used or the operator may use any two orthonormal 
unit vectors. Thirteen numbers for each vector 
must be inputted. If the sample is 10-dimensional, 
zeroes will be inserted for the last three 
components. A representation of the above example 
would look like: "А = 2919628, .31628,..., .31628, 
0,0,0" and a carriage return, "B = .31628,..., 
-.31628,0,0,0" and a carriage return, "*" and a 
carriage return. The XDS 9300 will compute the 
projections and relay the data to the AGI for 
display on the screen. The "X's" and "O's" denote 
the original "X's and Y's, respectively. 
(3) Interations and output 

(a) To change the orientation of the plane, the 
operator chooses either the "PICK X" or "PICK 
Y" function switches. The location of these 
buttons on the overlay are shown in Figure 
A-2-1. A video message, "TEXT EDIT REQUEST," 
will appear at the bottom of the screen. The 
operator depresses the button on the lightpen 
and places the light over the desired symbol 
(point) to be moved. Upon releasing the button, 
a portionof the letter wil] disappear to 


signify the selected point. The user pushes 
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(b) 


the two switches "COMPLETE CHOICE" and "END 
EDIT". The video message will disappear. The 
iteration will commence after one of the move- 
ment buttons (function switches) is depressed: 
ԱՆՆԱՆ UP "DOWN." The 9300 
recomputes the projections from the inner product 
with the orthonormal vectors and displays the new 
results on the screen. 

A printout of the new projections can be obtained 
by simultaneously depressing the "PRINT" and one 
of the movement function switches. The movement 
Switch must be released first. A copy of the new 
projections on the plane will be printed at the 
line printer. If the plot of the computer 
projections is desired, then depress the "PLOT" 


function switch. 
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APPENDIX 111. DATA LISTING 


This appendix presents a listing of the data sets used 
in the seven cases discussed in Section IV. The data has 
been arranged so that a row of numbers represents a single 
vector. For example, the data set for Case one is presented 
in three columns eens a three-dimensional entity. The 
thirteen columns in Case Seven represents a thirteen- 
dimensional problem. As an aid for reference purposes, the 


following index is provided: 


NUMBER NUMBER 


PAGE OF Y'S OF X'S 
CASE 1 21 ` 0 
CASE 2 20 20 
CASE 3 70 70 
CASE 4 40 20 
CASE 5 50 25 
CASE 6 38 44 
BASE 7 60 50 
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Բ о) го 


i o ^ ο са 


10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 


CASE 


ҮС 


140000 


1.0000 


2.5000 
2.5000 
4.0000 
4.0000 
140000 
1-0000 
895000 
2,5000 
4.0000 
4.0000 


1.0000 


140000 


4.0000 
4.0000 
2.5000 
197500 
1.7500 
3.2500 


3.2500 


1 
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DATA 


Ү(1,2) 


1,0000 
1.0000 
1.0000 
1.0000 
1.0000 
1.0000 
3.0000 
390000 
3.0000 
3.0000 
3.0000 
3.0000 
2,0000 
2,0000 
2,0000 
2.0000 
2.0000 
1.5000 
2,5000 
1.5000 


2,5000 


(173) 


1.0000 
2.0000 
1.0000 
2.0000 
1.0000 
2.0000 
1.0000 
2.0000 
1,0000 
2•0000 
1.0000 
2.0000 
1.0000 
2, 0000 
1.0000 
2.0000 
303229 
ce6614 
2.6614 
2.6614 


2.6614 





— هم ь. t‏ حم هم par e pa‏ هم 
о о Վ. OS лм о Гг +=‏ 


го 
O 


ο д о Վ ON O = о г + 


CASE 2 DATA 


Ү(1»1) 


390000 
4.0000 
5.0000 
340000 
4.0000 
5.0000 
5.0000 
3.0000 
3.0000 
4.0000 
5.0000 
390000 
390000 
4.0000 
5.0000 
5.0000 
5 «0000 
9.0009 
390000 


390000 


Y(I,2) 


2.0000 
2.0000 
2.0000 
4.0000 
4.0000 
4.0000 
3.0000 
3.0000 
2.0000 
2.0000 
2.0000 
3•0000 


440000 


4.0000 


4.0000 
3.0000 
2.0000 
4.0000 
4.0000 


2.0000 


Ү(1,3) 


5.0000 
5.0000 
5.0000 
50000 
5.0000 
5.0000 
5.0000 
5.0000 
7.0000 
740000 
7.0000 
7.0000 
7,0000 
7.0000 
7.0000 
7.0000 
6.0000 
6.0000 
6.0000 


6.0000 


10 
11 
ike 
13 
14 
13 
16 
17 
18 
19 
90 


CASE 
Х(ЈФ1) 


7,0000 
8.0000 
9.0000 
790000 
840000 
9.0000 
9.0000 
7.0000 
7 «0000 
840000 
9.0000 
7.0000 
740000 
8.0000 
9.0000 
9.0000 
240000 
9.0000 
7.0000 


7.0000 


లై 


65 


(CONT) 


Xx (INE) 


6.0000 
6.0000 
6» 0000 
8.0000 
8.0000 
8.0000 
740000 
7, 0000 
6.0000 
6. 0000 
6.0000 
7.0000 
8,0000 
8,0000 
8.0000 
7.0000 
60000 
8%0000 
80000 


6.0000 


X(Jo 3) 


9.000 
9.000 
9.000 
92000 
9.000 
9.000 
9.000 
9,000 
11,000 
11,000 
11,000 
11.000 
114000 
114000 
114000 
11.000 
10.000 
10.000 
10.000 


10.000 





4 м т 


Վ O տ 


10 
11 
12 
[e 
14 
15 
16 
17 
18 
19 
20 
ЕШ 


EC 


Ү(1»1) 


CASE 


Y(122) 


1.718 
1.621 
14599 
12-57 
1.328 
1.03% 
1.627 
1.452 
14255 
14120 
1.420 
14 365 
1.224 
1.996 
101291 
1.226 
12337 
14627 
1.648 
1% 175 
1:659 


J 532 


3 DATA 


Y(1,3) 


1.100 
1.526 
1.540 
15055 
1.639 
1.289 
14317 
14382 
1,711 
2.000 
14998 
162 095 
1.880 
1.541 
145278 
1.174 
1.319 
1.487 
1,962 
1615 
1643 


13937 


66 


YOTI) 


104892 
1,358 
14590 


1.176 


1,625 


1,646 
1812 
14452 
14215 
14072 
194861 
ass 


1.747 


14837 


1,638 
12806 
1,166 
1,687 
teres 
18577 
1.252 


1.162 


C145) 


1.617 
10552 
1.733 
1.892 
12754 
1.284 
1.197 
1.870 
1.546 
14344 
1.498 
1.460 
1.405 
1.046 
1,170 
1.244 
[1920 
1.422 
10378 
1.724 
1.842 


14 3/8 


Y(1,6) 


1,703 
1.746 
10484 
1,899 
1.684 
10520 
1.137 
10917 
14209 
1.397 
1.810 
1.825 
1,080 
1039 
1175 
1.802 
1.948 
1.413 
1.627 
1,216 
1,958 


1,827 


23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
44 
42 
43 


44 


Y(1,1) 


1,208 
1,968 
145490 
14692 
1,530 
19003 
14344 
1.649 
1.461 
1.701 
19745 
14386 
19757 
14122 
14289 
1.848 
14519 
1,312 


1,777 


“14504 


1.489 


1.599 


3 


Ел 


1.206 
1.342 
16153 
1,328 
1,909 
1.158 


1.059 


1.484 
1,181 
1,317 
1,655 
14039 
1,371 
11396 
156639 
1.17] 
1.047 
1.741 
19012 
1.072 
1,230 


1,176 


(CONT) 


Y(1,3) 


67 


ү(12,4) 


eco? 
1.487 
1౭౭9 
1.830 
12:22:77 
14096 
14204 
14244 
e397 
1,409 
14389 
14236 
186550 
1.647 
1.641 
1:775 
1.418 
2.528 
12232 
19236 
1593993 


19616 


Ү(1»5) 


1.497 
11523 
1,085 
1,106 
1,805 
1,783 
3515 
1.874 
1.924 
1.747 
1.368 
1.153 
14861 
1.294 
1,825 
1.545 
1,915 
1,055 
19442 
1.507 
1.287 


10422 


Y(1,6) 





CASE 3 (CONT) 


45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 


68 


(171) (176) Y(1,3) У(12%) Ү(1»9) Y (126) 
14997 14119 14400 10364 14144 1,194 
1,239 1,355 14864 14233 14052 1.657 
1,842 14195 14941 1176 1.052 1.950 
ποι 14840 1.602 sas 1.907 1,756 
14007 1.528 1.540 1.648 1,532 14372 
14952 1.902 1.775 14218 1.930 1.344 
1.668 1.392 10636 ee 13292 1.987 
esse 1.263 1.838 14016 14905 1.538 
1.492 1.989 14514 12500 1.607 1.803 
1.629 14333 14027 14080 1.741 1.282 
14977 1.656 1.145 1.948 1372 {2555 
14474 1.988 1.088 1.578 1.882 1.557 
1.904 1.334 14104 1.974 1.088 1.021 
1.184 1.770 1.259 terre 1.551 = 
1.029 ieee 1.112 14501 12411 1,352 
14582 1,701 14598 СБА 1.025 14647 
1.592 1.496 =e 14020 1.805 12765 
1.296 1.101 WEIT 14092 1.789 1,297 
1,059 1.146 1.531 1.058 1.992 1.649 
14165 1,174 1.485 10724 14140 1.005 
1.430 1,907 1,089 14360 14180 1.028 
1.746 1.939 1.344 ΠΕ 14626 12260 





67 
68 
69 
70 


ΥΓ 


1.203 
1.431 
1.286 


1.089 


3 


Noe 


1.601 
1.050 
1.170 
1.497 


(CONT) 
Y(1,3) 


1.525 
1.242 
1.939 


1.623 


Ү(1,%) 


14610 
1.270 
19011 
14272 


Y(1,5) 


1,575 
1.146 
1.486 


14532 


Y(1,6) 


14151 
1.874 
1.768 


1.876 





ա గు 


ч ο (щл Հ 


10 
11 
ΤΕ 
13 
14 
15 
16 
17 
18 
19 
20 
21 
ԷՔ 


Х(Ј,1) 


«528 
«063 
#2541 
«466 
293 
«023 
239 
4402 
«335 
«573 
9271 
4949 
967 
4928 
‚112 
694 
‚171 
0543 
233 
4839 
4182 


331 


CASE 3 (CONT) 


X(Js2) 


‚236 
843 
785 
‚538 
«315 
‚178 
‚141 
‚277 
‚318 
‚189 
«905 
.428 
+284 
• 429 
(149 
«080 
‚265 
«909 
«987 
«327 
.594 


«007 


Х(.93) 


e840 


799 


‚924: 


4899 
990 
• 375 
505 
699 
4463 
2085 
«657 
0674 


• 305 


e243 


4884 
420 
«374 
‚612 
43%6 
‘231 
(759 


«058 


70 


Х (125) 


«390 
4273 
4509 
4784 
«536 
0947 
4949 
66% 
«551 
«751 
‚175 


.400 


.450 


.477 


582 
715 
‚528 
064 
• 305 
ος 
«759 
e 455 


Х(, 5) 


4761 
480% 
. 571 
«Jol 
‘209 
«090 
«500 
«139 
«079 
• 155 
‚275 
4800 
688 
703 
MESE 
«068 
e462 
9576 
748 
. 198 
e420 


«072 


X(Js6) 


4516 
• 027 


571 





J 


23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 


44 


X( Jel) 


• 366 
«067 
. 304 


1.000 


270% 


4702 


` «506 


•766 
•038 
4683 
296 
|, 337 
• 398 
«476 
0477 
«166 
«596 


366 


“4060 


4833 
6355 


4112 


CASE 


3 (CONT) 
X(Jo 2) Χ(ώ, 3) 
942 747 
288 • 326 
4483 «619 
‚257 «693 
• 505 “329 
4388 «307 
«967 «698 
օ787 458% 
«807 «826 
0534 «796 
ebi2 «109 
«676 «070 
«659 29115 
«091 «688 
392 «514 
«004 «062 
41102 «665 
‚208 «867 
4130 e542 
«373 «658 
9541 e412 
436% ‚799 


7] 


Х(Ј»9) 


010 
4129 
«135 
«116 
«180 


«345 


‚029 
131 
0044 
098 
.692 
,882 
‚768 
‚552 
.633 
‚839 
0498 
(493 
«230 
0486 
‚861 


«807 


X(J55) 


(371 
0143 
4713 
419% 
«806 
‚351 
«102 
«092 
«988 
‚963 
‚836 
0448 
‚378 
‚897 
«048 
(412 
‚376 
‚592 
«975 
‚85% 
‚513 


• 380 


Х (416) 


(237 
(353 
.821 
,278 
718 
eS 
‚656 
5119 
‚259 
‚982 
• 105 
305 
1627 
,989 
520 
591 
Soul 
135 
5796 
266 
,525 


“097 





B 


45 
4 6 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 


66 


X( Jal) 


e325 
413% 
9767 
234 
«318 
«890 
«621 
4521 
«718 
826 
940 
«636 
4780 
«984 
.884 
AT 
«165 
«558 
«570 
2511 
559 


«^03 


CASE 


3 (CONT) 
Χ(ν, 8) Χ(ω, 3) 
‚321 ‚120 
692 «961 
734 • 132 
«679 • 345 
• 306 ‚583 
4403 (295 
5712 0449 
105 0245 
1146 622 
«635 266 
«189 . 33 
0978 629 
‚202 «703 
(188 «637 
289 «694 
P" 901 
0094 0425 
0528 «688 
‚175 4717 
066 4264 
(044 41/5 
931 299 


72 


«(աձ) 


010 
4127 
436 
«683 
0455 
‚126 
. 158 
«684 
‚702 
‚808 
‚178 
• 380 
«07e 
«680 
«3^0 
410% 
0477 
nos 
6 
498 
«705 


.654 


X(J,5) 


• 338 
•770 
0145 
4080 


4137 


2138 _ 


«108 
«334 
9523 
«058 
9540 
6292 
4150 
»416 
.876 
‚638 
698 
953 
• 320 
872 
«016 


«600 


X( Jo» 6) 


4971 
188 
476 
9350 
«658 
1222 
‚278 
‚728 
ΖΡ 
4993 
4829 
4874 
«824 
«818 
«109 
«232 
4048 
e264 
‚253 
‚186 
«562 


e 733 





‚т 


67 
68 
69 


70 


Χ( 1) 


6463 
«739 
e50^ 


‚212 


СА$Е 


3 (CONT) 
«(ա2) ౫ (౮౫౩) 
.481 ‚116 
«094 0651 
| «663 «638 
| «056 ‚592 


73 


(աձ) 


.Օ7Կ 
.967 
‚712 


‚658 


X(J,!5) 


«870 
9349 
5513 


.391 


Х(,96) 


90 1 3 
e173 
822 


«566 





c! + QQ M ౯ 


о 


10 
11 
Je 
13 
1^ 
15 
16 
17 
18 
13 
20 
21 
се 


САЗЕ 


ү(1,1) 


"0500 
“1,2000 
ге 5000 
1,0000 
1,0000 
1,2000 
-• 5000 
“90590 
1,1000 
1,1000 
"1.1000 
„1,1000 
1,5000 
1,5000 
©0470 
“1,2000 
”1,2000 
1,0000 
1,0000 
e2000 
=». 682C 


•• 7650 


4 


DATA 


Ү(122) 


142000 
2.0500 
«1.0000 


". 7000 


.5000 


«7000 
140000 
1.0540 
1.1000 

"1.1000 
1.1000 
-1+.1000 

.5000 

-?5000 
„1,5000 

«5000 
»е 5000 
4000 
”4000 
«51000 
63596 


1.4370 


74 


ΥΓ» 


«7000 
-1.2000 
-.5000 
+6000 
1.0000 
- 7000 
«6000 
1,1500 
«6000 

e ¿2000 
«6000 
-,2000 
«8000 

• 7500 

= |+ 3000 
1.0000 
140000 
- 1•%000 
- 1.4000 
1,0000 
• • 1320 


‚6250 


ХМ (124) 


10000 
40500 
-1.0000 
1.0000 
«8000 
-1.1000 
1.2000 
«0550 

• 2000 

• 2000 
"0000 
-.2000 
«6000 
7500 
«9400 
«5000 
"5000 
«4000 
..4000 
1.3000 
“1,8380 


«8080 





23 
e^ 
25 
26 
27 
28 
29 
30 
31 


32 


33 
34 
35 
36 
37 
38 
39 


40 


CASE 
Ү(1,1) 


«5660 
#¢7580 
"1,9940 
«7060 

• 3690 
”.84%0 
» 7320 
"7350 
“1,9130 
«8180 
=+ 3800 
% 5510 
“1,9410 
41130 
0980 
“1,9060 
”190050 


„1,9770 


4 


(CONT) 
Y(1,2) 


•. 4680 


"25810. 


0860 
#07940 
-.8250 
#02490 

«0710 

«0330 

©3420 
1.4470 
1.4660 
1.5800 
3660 
1,5680 
1,5600 

0180 
«4050 


- 0900 


75 


(173) 


= + 2980 
1.7140 
•. 1110 
1,6710 
1,6220 
«0820 
«53260 
#20700 
«1820 
“,9770 
“1,1350 
“+9870 
• 1630 
«1.0680 
«3780 
-. 4960 
1.6320 


.0450 


Ү(199) 


«128480 
.3730 
«0550 
• 2390 
«5890 
"1,7860 
"197740 
"18580 
0630 
©0940 

- ©6260 
«4580 
«0220 
«3450 
-1.1680 
«1410 
41690 


-.0510 





00 ՀՎ о ся + а) го -- 


O A 


12 
13 
14 
15 
16 
17 
18 
19 


20 


CASE 


X(Ja1) 


«1370 


2450 
1140 


48360 
«1970 
«0770 
«6170 
«0970 
«2930 


4 


«6080 
«1990 
«2670 
«0680 


44860 


“0540 
41630 
«4100 
42020 


.3550 


3 


«6580 


5 1 


4 


(CONT) 


X(Js2) 


«6680 
7640 
г«8740 

«5580 

«5910 
- «9320 

«6550 

“S710 

*6520 
«7060 
- «8330 
«7540 
• 9320 
"7850 

«6090 
= 8880 

«3390 
7230 

«5350 


«5960 


76 


X(J,3) 


= 5700 
"5690 
, 2970 
4460 
«4740 
"2750 
"6260 
” (7950 
©2490 
~ «4710 
«3580 
• 3010 
- «3580 
«1990 
«2680 
e2430 
(3970 
5960 
- 6690 


. 4730 


24274) 


e e4500 
-+0990 
• 3620 
«2870 
5380 
«1710 
«1710 
«1160 
«2060 
.օ3580 
21330 
«2510 
= 0880 
«2100 
-.5980 
+2220 
«3460 
«1500 
«3320 


-.5370 





ғ о го 


Վ ԺԹ Մ 


10 
11 
12 
18 
14 
15 
16 
17 
18 
1- 
20 
21 
22 


CASE 


Ү(1,1) 


5.1000 
4.9000 
4.7000 
4.8000 
5.0000 
5.4000 
4.6000 
5.0000 
4.8000 
4.9000 
5.4000 
4.8000 
4.8000 
4.3000 
5.8000 
5.7000 
5.4000 
5,1000 
5.7000 
5.1000 
5.4000 


5. 1000 


5 DATA 


212112 


395000 
3.0000 
3,2000 


3.1000 


3%6000 


3.9000 
3.4000 
3.4000 
2౪9000 
3.1000 
3.7000 
3.4000 
340000 
340000 
4.0000 
4.4000 
3е9000 
3.5000 
3.8000 
3.8000 
3.4000 


3.7000 


77 


У Ки З) 


1.4000 
1.4000 
1,3000 
1,5000 
1,4000 
1.7000 
1.4000 
145000 
1,4000 
145000 


1,5000 


1.6000 


1.4000 
1:1000 
142000 
1.5000 
1» 3000 
1.4000 
1.7000 
1,5000 
1.7000 


195009 


Ү(1,%) 


42000 
©2000 
«2000 
©2000 
«2000 
©4000 
«3000 
.2000 
22000 
«1000 
«e000 
©2000 
«1009 
• 1000 
•2000 
©4000 
«4000 
«3000 
«3000 
«3000 
“2000 


«4000 





23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 


4 


CASE 5 


Ү(191) 


4.6000 
5.1000 
4.8000 
5.9000 
6,1000 
6,3000 
6.1000 
6•%000 
7«0000 
6•%000 
6•9000 
5.5000 
6.5000 
5.7000 
6• 3000 
4.9009 
6.6000 
5,2000 
5,0000 
5.9000 
6.0000 


6: 1000 


(CONT) 


Y(1,2) 


3.6000 
3.3000 
3.4000 
3.2000 
2.8000 
2.5000 
2.8000 
2.9000 
3,2000 
3,2000 
3,1000 
2,3000 
2,8000 
2.8000 
3.3000 
2.4000 
893000 
6.7000 
0.0000 
3.0000 
2.2000 


893000 


78 


ҮСРЭЭ! 


1.0000 
1.7000 
1.9000 
4.8000 
4.0000 
4.9000 
4.7000 
4.3000 
4.7000 
4.5000 
4.9000 
4.0000 
4.6000 
4.5000 
4.7000 
3.3000 
4.6000 
349000 
3.5000 
4.2000 
4.0000 


4.7000 


Ү(1,%) 


©2000 

«5000 

e2000 
1.8000 
1,3000 
1,5000 
1,2000 
1,3000 
1.4000 
1.5000 
1.5000 
193000 


145000 


_ 1.3000 


146000 


1.0000 


163000. 


1.4000 
1.0000 
1.5000 
1.0000 


144000 





45 
46 
47 
48 
49 
50 


CASE 5 


Ү(1,1) 


5.6000 
6.7000 
5.6000 
5.8000 
6.2000 


5.6000 


(CONT) 


(176) 


2.9000 
3.1000 
3.0000 
2.7000 
042000 


045000 


γην) 


3.6000 
4.4000 
4.5000 
4.1000 
4.5000 


3.9000 


Ү(1 28) 


1.3000 
144000 
195000 
1.0000 
1.5000 


141000 





10 
11 
де 
13 
15 
15 
16 
17 
18 
19 
20 
el 
Be 


CASE 5 


X( Jal) 


6.5000 
6«4000 
6«8000 
5.7000 
5.8000 
6.4000 
6.5000 
7.7000 
7.7000 
6.0000 
6. 3000 
5.8000 
7,1000 
6• 3000 
6.5009 
7.6000 
4.9000 
7.3000 
6.7000 
7.2000 
6.9000 


5 «6000 


(CONT) 


Х(, 2) 


3.2000 
2,7000 
390000 
2.5000 


2,8000 


3.2000 
3.0000 
3,8000 
2.6000 
2,2000 
393000 
2,7000 
390000 
2.9000 
390000 
3«0000 
25000 
2.9000 
2.5000 
3,6000 
3.2000 


2.8000 


80 


5) 


5.1000 
53000 
5.5000 
5.0000 
5.1000 
5.3000 
5.5000 
6.7000 
6.9000 
5.0000 
640000 
5,1000 


5,9000 


5.6000 


5.8000 
626000 
4.5000 
6. 3000 
5.5000 
6.1000 
5.7000 


4.9000 


X(Ja ^) 


2.0000 
199000 
2౪1000 
ге 0000 
2.4000 
2.3000 
1,8000 
2,2000 
2,3000 
145000 
2.5000 
1.9000 
241000 
18000 
2.2000 
2,1000 
1,7000 
1,8000 
1,8000 
2.5000 
2.3000 


2۰0000 





23 
24 
po 


BASE 5 


Χ( 1) 


7.7000 
6.3009 
67000 


(CONT) 


Xx(Ja2) 


2.8000 


997000 
393000 


8] 


Х(. 3) 


6 7000 
& 9000 


5.7000 


XtJo ^) 


20000 
1.8000 
2.1000 





I 


ω oO ο 


о AO O Վ ο οι Հ 


Te 
13 
14 
15 
16 
17 
18 
19 
20 
21 


ве 


Y(1,1) 


160. 
155 
130. 
161, 
135. 
165. 
150» 
1480 
150. 
120. 
1206 
190» 
100, 
150, 

90» 
120. 
200» 
190. 
105. 
C106 

90» 


1120. 


(176) 


6 DATA 


CASE 
YCI,3) Y(I,4) 
10.0 08.0 
8.0 27,0 
1040 26.0 
10.0 67.0. 
12.0 27.0 
11.0 23,0 
9.0 29.0 
8.0 26,0 
7.0 3140 
6.0 33.0 
8.0 32.0 
8.0 3090 
9.0 35.0 
7«0 69.0 
8.0 3040 
7.0 35.0 
9.0 30«0 
8«0 28.0 
7.0 27.0 
8.0 26.0 
10.0 25,0 
11160 25۰0 


82 


Y(1,5) Y(1,6) 


б 
85, 
75, 
70. 
88. 
95. 
90. 
85. 
60. 
BE 

m 
55, 


480 


65. ° 


70. 
67. 
62, 
63. 
65, 
67. 
68. 


60», 


665. 
570. 
675. 
580. 
390. 
435. 
440. 
650» 
640, 
4306 
650. 
655. 
645. 
660. 
530. 
435. 
440. 
430 


530% 


Ү(1»7) 


Ү(1»8) 


61.0 
17.0 
19.0 
61.0 
03.0 
28.0 
25.40 
28.0 
60.0 
18.0 
23.0 
27.0 
26.0 
20,0 
17*0 
60.0 
63«0 
19.0 
27.0 
2040 
2140 


22.0 





I 


23 


25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 


38 


100» 
70. 
100. 
130. 
90. 
80. 
Обе 
70. 
80. 
95. 
100. 
95, 
85, 
70. 
80 є 


55. 


443 
4.5 
4.8 
ове 


4.5 


е 


449 


442 
541 
446 
448 
це? 
408 
5.4 


5.5 


CASE 6 


Ү(1,1) Ү(192) Ү(1,3) 


9.0 

8.0 

9.0 

9.0 
11.0 
1040 
10,0 
1ге0 
1340 
1570 
11.9 
10.0 
12.0 
11,0 
10,0 
13,0 


2590 
23,0 
2790 
25,0 


37,0 


32,0 


25,0 
3040 
36,0 
3140 
24.0 
2740 
65.0 
Շ6.0 
2140 


3340 
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(CONT) 


70. 
6%. 
65. 
0% 
744 
78， 
82. 
85, 
69» 
70» 
76. 
7%. 
73. 
78. 
80. 


81. 


ТҮСІМІ? ҮСІ»5) ҮСІ;6) 


8840. 
4500 
355. 
380. 
350. 
450. 
260. 
262, 
265, 
265. 
2070. 
955 


360» 


365. 


370. 


252% 


МАУ) 


(178) 


1740 
19.0 
2440 
29.0 
31.0 
32,0 
23.0 
2440 
66.0 
23.0 
39.0 
31.0 
69.0 
3040 
30.0 


35.0 





J 


O (Q + оу го 


О O ы 


11 
te 
13 
14 
15 
16 
17 
18 
19 
20 
21 
ec 


200. 
260. 
195. 
195. 
220. 
880. 
190. 
285. 
300. 
225. 
260. 
280» 
300. 
310. 
290. 
260» 
290. 
160. 
240. 
195. 
290. 


610. 


4(1) 214126) 


304 
3.1 
3.0 
3.2 
3.3 
3.0 
3%% 
3.0 
3.0 
3.0 
304 
390 
Зей 


3.0 


3.0 
3.3 
3.1 
3.5 
3.1 


304 


CASE 


US Y 


10.0 
8.0 
9.0 
990 

10.0 


9.0 
11.0 
9.0 
10.0 
8.0 
8.0 
10.0 
1149 
1240 


9.0 
1190 
1140 

8.0 
10.0 


9.0 


6 


2440 
21.0 
20.0 
1940 
2440 
25.0 
26.0 
19.0 
20,0 
22۰0 
22۰0 
20۰0 
20۰0 
19.0 
26.0 
22.0 
౭50 
20.0 
2040 
2140 
1940 


22,0 


84 


(CONT) 


98. 
110. 
105. 
110. 

95, 

90» 

96. 
100. 
108. 
108. 

97. 
112. 
108. 
106. 

94. 

98, 
100. 

79. 

88. 
En 

TA 


96, 


X(J,#&) xtJa5) X(Ja 6) 


210. 
220° 
130. 
010. 
205» 
810. 
070. 
990. 
120. 
985. 
090. 


200. 


835. 


055: 
240% 
015. 
010. 
170. 
990. 
975. 
860. 


950. 


X(Js7) 


X(Js8) 


40.0 
360 
4400 
36.0 
29.0 
34.0 
30.0 
31.0 
4090 
25.0 
38.0 
41.0 
45.0 
37.0 
8140 
29.0 
32.0 
34 + 0 
42.0 
44.0 
40.0 


47.0 





J 


23 
24 
2 
26 
27 
28 
29 
30 
31 
зе 
H3 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 


4 


1806 
205. 
ԲԱՑ: 
270», 
290, 
320, 
210. 
210. 
185. 
200», 
170. 
140. 
110» 
100» 
95. 
65. 
55 
70. 
85 є 
115: 
IR 


1200 


mig XU 2 ) 


CASE 6 


22.0 
23,0 
61.0 
20,0 
03.0 
19.0 
18.0 
21,0 
25,0 
86.0 
60.0 
60.0 
22,0 
1940 
2440 
3040 
07.0 
08.0 
19.0 
2140 
67.0 


24.0 


85 


(CONT) 


97. 
90. 
100. 
111。 
106. 
87. 
116. 
99 
96: 
98. 
95e 
98: 
130. 
1280 
164. 
134, 
128• 
118, 
117. 
122, 
114. 


120. 


321 Խ4)1721125) "1176) 


990» 
805. 
700. 
170. 
350. 
160. 
010% 
190» 
055. 
980. 
095. 
990, 
2020 


216. 


218. 


200» 
205. 
208, 
206. 
198. 
ΘΟ ὃ 


24%. 


Χίν 7) 


X(Js8) 


28.0 
38.0 
32.0 
31.0 
33.0 
38.0 
45.0 
8040 
41.0 
38.0 
40.0 
32.0 
19.0 
1460 
10.0 
12.0 
11.0 
12.0 


1440 





ғ оо గు 


ο о 


ՀՎ 


11 
T2 
13 
14 
15 
16 
17 
18 
i9 
20 
21 
22 


Y1 


Ye 


365 
τ 
343 
432 


376 


410 
408 
410 
315 
546 
437 
506 
338 
446 
426 
410 
483 
377 


410 


363 


335 


τ» 


152 
126 
180 
110 
115 
110 
105 
150 
118 
ο. 
120 
106 
110 
120 
135 
140 
120 


130 
1720 
140 


186 


Y^ 


20 
70 
86 
65 
80 
70 
65 
100 
76 
70 
70 
Z 
70 
70 
80 
85 
70 
70 
85 
80 
90 


82 


CASE 7 DATA 


Y5 yé 
4 2 
1 2 
1 1 
ES 
216! 
о 
1 1 
1 2 
1 2 
4 1 
о 
4 2 
2 2 
1 4 
1 1 
TOM 
1 e 
2 2 
2 і 
4 2 
. 1 
3 1 


86 


Υ7 


Q ա. Թ 


Y8 


то ә го го о о 


го о го 


гә го 


то ә ә , гә го 


Y3 


го m 


то vv m 


MO Үү Ге 
1 2 290 
2 2 356 
3 С» 4276 
4 с 2278 
3 Өз “276 
3 ౭.6౭55 
3 со 276 
4 атаа 
& 593.5 
3 12222182 
s 2 3/2 
3 17221533 
З 199278 
4 17222278 
4 1 1275 
ч 1 250 
e 1. 279 
4 5 2390 
4 253269 
3 5 246 
4 11278 
© 12278 


ὙΠ 


39 
38 
42 
42 
41 
43 
50 
44 
28 
47 
39 
37 
39 
50 
48 
35 
48 
38 
39 
49 
38 


38 





23 
24 
242 
26 
27 
28 
mo 
30 
31 
зе 
ՅՅ 
34 
85 
36 
37 
38 
39 
40 
41 
42 


43 


44 


Y1 


Ye 


410 
422 
517 
250 
466 
366 
271 
467 
372 
300 
348 
288 
477 
406 
489 
406 
367 
444 
402 
462 
ST 


410 


Y3 


140 
152 
110 
130 


120 
110 
130 
128 
140 
170 
130 
115 
140 
110 
130 
195 
125 
122 
120 
140 


140 


Y4 


80 
Эс 
80 
70 
74 
90 
70 
70 
Эс 
75 
80 
80 
80 
110 
64 
80 
85 
80 
70 
78 
25 


90 


BASE 7 


Y5 


h 


Y6 


87 


(C8NT) 
Y7 Y8 
3 1 
e e 
3.072 
3 2 
e e 
e e 
3 e 
3 2 
3 e 
3 e 
2 1 
& 2 
2 2 
27 
ο ο 
4 2 
3 2 
3 2 
1 e 
3 1 
e e 
e e 


Y3 


гу 


го ә ә ә гә т го 


го 


Խ г гә гу 


ածն 111115 
4 e 268 
3 24-278 
4 э 250 
4 6.5 
4 3 ps 
4 3,229 
4 б. 2295 
4 9, 2278 
ч D XS 
4 5 200 
4 221335 
ч e 268 
4 2 25 
4 5 300 
З 5 278 
3 ӨН 2275 
3 3 ο 
з 1 зог 
3 tees 
4 1 278 
З 1 350 
3 e 245 


Y13 


40 
49 
53 
45 
38 
47 
36 
41 
43 
38 
40 
47 
43 
33 
42 
51 
40 
45 
46 
41 
45 


Эс 





45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 


Yi 


Ve 


446 
468 
410 
410 
355 
394 
312 
383 
489 
410 
404 
378 
369 
388 
483 


461 


ҮЗ 


уч 


60 
70 
80 
88 
70 
80 
60 
60 
80 
82 
98 
90 


70 


90 
80 
90 


CASE 7 


Y5 Y6 
1 1 
1 e 
1 2 
1 2 
1 1 
1 2 
1 2 
1 2 
1 2 
1 E 
1 2 
1 e 
e 1 
3 1 
4 1 
1 6 


88 


(CONT) 
GE 
4 1 
3 2 
3 с 
3 e 
3 e 
133117 
1 
1 6 
e 2 
г e 
e 8 
3 1 
4 2 
S 
T 
4 1 


Y9 


го ә ә гә € 


го 


ΟΥΤΕ 
3 e 202 
4 17 313 
5 5. 351 
4 1 242 
4 3. 366 
4 е 285 
4 e 310 
4 2 310 
ఓ 2 278 
L 0: 6-0 
& e 278 
3 1 267 
3 1 278 
4 5 288 
4 27] 
e 14226 


113 


40 
37 
40 
це 
23 
42 
47 
47 
46 
41 
43 
43 
44 
43 
39 
52 





ν г 


Xe 


404 
418 
288 
349 


396 


411 
443 
479 
388 
445 
389 
406 
431 
413 
411 
406 
444 
375 
383 
439 
489 


387 


x3 


160 
118 
ites 
170 
118 


180 


X4 


78 
76 


86 


110 


7e 


120 


90 
90 
90 
88 
78 
70 
60 
86 
50 
76 
80 
70 
80 
80 
70 


90 


CASE 7 (CONT) 


mm о го 
о ә  σ гә го 


ον го го го 


X10 X11 x12 


209 
298 
209 
343 
сеї 
196 
271 
260 
681 
215 
213 
248 
се: 


247 


248 
151 
248 
236 
248 
221 


273 


67 
43 
37 
40 
37 
Se 
37 
40 
40 
42 
32 
40 
34 
41 
40 
23 
42 
41 
36 
40 
4 4 


43 





23 
24 
25 
26 

27 
ET 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 


42 


43 


44 


х1 


х2 


351 
350 
536 
432 
393 
369 
418 
383 


325 


395 


375 
432 
267 
276 
489 
369 
398 
442 
394 
444 
389 
359 


ХЗ 


120 
110 
124 
110 
150 
132 
120 
126 
155 
120 
117 
142 
120 
158 
188 
150 
150 
118 
130 


120 


130 


x4 


90 
66 
84 
70 
90 
80 
80 
84 
80 
78 
75 
80 
84 
100 
70 
84 
90 
70 
80 
78 
95 


100 


CASE 7 


X5 X6 
4 2 
2071 
. 2 
TNCS 
ՅՆԱ 
1 2 
3 2 
14 ὁ 
1 e 
E 
125! 
. 2 
4 2 
6 
TE 
2071 
. 2 
3 e 
1 e 
3 e 
ЁС 
го 


90 


(CONT) 
х7 X8 
3 3 
3 3 
2 3 
4 1 
з. 
2 3 
3 3 
4 3 
4 3 
H 1 
4 

4 3 
4 

4 

3^ 

4 3 
4 3 
e 3 
3 3 
e 3 
° 3 
2 3 


ա Թա. о 


29 


X10 X11 x12 


258 


248 
282 
240 
285 
227 
236 
221 
175 
178 
ος 
233 
213 
167 
248 


294 


223 


240 
288 
256 
273 
155 


X13 


36 
41 
33 
55 
46 
40 
40 
36 
30 
40 
39 
43 
40 
41 
44 
43 
42 
44 
31 
42 
ць 


35 





45 
46 
47 
48 
49 
50 


х2 


ese 
379 
378 
381 
304 
437 


«3 


142 
106 


140 


199 
128 
140 


X^ 


100 
72 
100 
20 
82 
90 


X10 x11 x12 


ез 
264 
309 
све 
250 
248 


x13 


55 
34 
39 
54 
33 


36 





APPENDIX IV. 


TABLE OF A AND B VECTORS 


The following table is the values of the various components 


of the A and B vectors used in the seven examples. 


of the thirteeen components have been listed. 


To use the 


Only six 


table, repeat the pattern of the values until the number of 


components equals the dimension of the problem. 


For the 


remaining components, use a zero value until one has filled 


the thirteen components. 


DIMENSION 


11 


12 


(2 > 00 > Ջա 2 ա 2 Ὁ > ա > (9 > (9 > (2 2> (2 > wo m 


TABLE OF A AND B VALUES 


. 50000 
.50000 
. 90000 
. 50000 
.40825 
‚40825 
‚40825 
‚40825 
295595 
239355 
229059 
35355 
21028 
201225 
31628 
. 31628 
.28867 
‚288657 
.28867 
.28867 
„26/26 
‚26726 


‚70711 
ΠΠ] 
. 50000 
.50000 
. 50000 
. 50000 
. 40825 
. 40825 
. 40825 
. 40825 
. 35555 
ο 
19 9055 
230399 
31028 
„1628 
31628 
“31528 
223667 
229907 
.2886/ 
.28867 


COMPONENTS 
3 4 
.50000 .00000 
.50000 .00000 
.50000 .50000 
‚50000 -.50000 
‚40825 .50000 
.40825 -.50000 
‚40825 .40825 
‚40825 -.40825 
‚35355  .40825 
‚35355 -.40825 
‚35355 .35355 
‚35355 -.35355 
‚31628  .35355 
Լւ 4368355 
.31628 «31628 
‚31628 -.31628 
.28867 .31628 
.28867 -.31628 
.28867 .28867 
.28867 -.28867 
,26726 .28867 
.26726  -.28867 


.00000 
.00000 
.00000 
.00000 
‚40825 
‚40825 
‚40825 
‚40825 
225299 
299315 
239999 
35359 
31028 
2931629 
31028 
21625 
-29867 
.28867 
728207 
222881 
. 26726 
.26726 


.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.40825 
.40825 
‚40825 
‚40825 
225355 
399350 
215315 
III 9 
626 
231623 
Эле 
221525 
"28867 
. 28867 
225552 
noob 





APPENDIX V. TWO DIMENSIONAL PLOTS 


- The following figures represent the plots of the computer 
printout. The values of the A and B vectors are printed below 
the plot if a discriminant function was determined. The 
components for the initial projections can be extracted from 
Appendix IV. A polaroid print is also included for each case. 


This allows for a comparison between the AGT display and the 


computer plot to be obtained. 
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Figure IV-1-1, CASE 1, STEP 1 
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Figure IV-l-2. CASE, IT, STEP2 
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Figure IV-1-3. Case 1, STEP 3 
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Figure IV-1-4. CASE 1, STEP 4 
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Figure IV-3-4. CASE 3, STEP 4 
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Figure IV-3-5. CASE 3, STEP 5 
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Figure IV-4-1. CASE 4, RUN 1, STEP 1 
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Figure IV-4-2. CASE 4, RUN 1, STEP 2 
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Figure IV-4-3. CASE 4, RUN 1, STEP 4 
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Figure IV-4-4. CASE 4, RUN 1, STEP 5 


116 








Figure IV-4-4P. CASE 4, RUN 1, STEP 5 
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Figure IV-4-5. CASE 4, RUN 1, STEP 6 
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Figure IV-4-6. CASE 4, RUN 2, STEP 1 





ΑΦ ΕΣΑΚΕ ΝΑΣΑ ΝΑ ΚΑΛΑ ΝΑΟΣ ΚΑ ΟΣΕΡ HE OE HE аж 


Ж Ф 
: # 
4 з 
5 4 
: * 
: 4 
. 0 E 
x # 
* # 
- 0 * 
« Օ 0 + 
* " * 
4 Х 0 0 5 
4 ЕХ E O ο 0 5 
4 0 O = 
- X ο ο * 
. Ж ОХ κοπο ος ο 0 . 
4 Ё ХЕ 0 Ж 
» 0 0 * 
# 0 * 
E 4 
# 0 5 
8 000 τ 
4 0 D O ο е 
Š 4 
5 4 
L 4 
` * 
E * 
E శ 
8 4 
* 0 . 
8 00.0 * 
+ 20 % 
қ ЕКЕ кВт кекете + 


Figure IV-4-7. CASE 4, RUN 2, STEP 2 
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A(3) = .3773 A(7) = .3191 B(3) =-.0270 B(7) = .1666 
A(4) = .2770 A(8) = .3455 B(4) = .3068 B(8) = .0788 
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Figure IV-4-8. 
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Figure IV-4-9. CASE 4, RUN 2, STEP 4 
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Figure IV-4-10. CASE 4, RUN 2, STEP 6 
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Figure IV-4-10P. CASE 4, RUN 2, STEP 6 
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Figure 1V 5-1. CASE 5, STEP 1 
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Figure IV-5-2. CASE 5, STEP 2 
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Figure IV-5-5. CASE 5, STEP 6 
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Figure IV-6-2. CASE 6, RUN 1, STEP 2 
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Figure IV-6-3. CASE 6, RUN 2, STEP 4 
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Figure IV-6-4. CASE 6, RUN 2, STEP 1 
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Figure IV-6-4P. CASE 6, RUN 2, STEP 1 
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Figure IV-7-2. CASE 7, RUN 1, STEP 2 
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Figure IV-7-3. CASE 7, RUN 1, STEP 11 
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Figure IV-7-4. CASE 7, RUN 1, STEP 15 
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Figure IV-7-5. CASE 7, RUN 1, STEP 20 
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Figure IY-7-7. CASE 7, RUN 2, STEP 1 
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Figure IV-7-8. CASE 7, RUN 2, STEP 2 
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Figure IV-7-9. CASE 7, RUN 2, STEP 9 
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Figure IV-7-10. CASE 7, RUN 2, STEP 12 
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Figure IV-7-11. CASE 7, RUN 2, STEP 14 
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Figure IV-7-12. CASE 7, RUN 2, STEP 27 
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